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Abstract. In this study we discuss two key issues related to a small-scale dynamo 
instability at low magnetic Prandtl numbers and large magnetic Reynolds numbers, 
namely: (i) the scaling for the growth rate of small-scale dynamo instability in the 
vicinity of the dynamo threshold; (ii) the existence of the Golitsyn spectrum of 
magnetic fluctuations in small-scale dynamos. There are two different asymptotics for 
the small-scale dynamo growth rate: in the vicinity of the threshold of the excitation 
of the small-scale dynamo instability, A oc ln(Rm/Rm CI ), and when the magnetic 
Reynolds number is much larger than the threshold of the excitation of the small- 
scale dynamo instability, A oc Rm 1 ' 2 , where Rm cr is the small-scale dynamo instability 
threshold in the magnetic Reynolds number Rm. We demonstrated that the existence 
of the Golitsyn spectrum of magnetic fluctuations requires a finite correlation time of 
the random velocity field. On the other hand, the influence of the Golitsyn spectrum 
on the small-scale dynamo instability is minor. This is the reason why it is so difficult 
to observe this spectrum in direct numerical simulations for the small-scale dynamo 
with low magnetic Prandtl numbers. 
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1. Introduction 

Generation of magnetic field by turbulent motions of conducting fluid is a fundamental 
mechanism of magnetic fields observed in stars, galaxies and planets. There are different 
kinds of turbulent dynamos: large-scale and small-scale dynamos. The large-scale mean- 
field dynamo implies that the amplification of magnetic field occurs at scales which are 
much larger than the maximum scale of the turbulent motion. This kind of dynamo 
includes: (i) the aQ and a 2 fl dynamos caused by the combined action of the a effect 
and differential rotation (see, e.g., [TJ [2J [3J, SJ [5]); (ii) a 2 dynamo in helical turbulence; 
and (iii) the shear dynamos in non- helical turbulence [6j [TJ E] . 

On the other hand, generation of magnetic fluctuations occurs at scales which 
are smaller than the maximum scale of the turbulent motions (see, e.g., reviews 
0, Uni EEIJ 121 EEl EH]). Self-excitation of magnetic fluctuations with a zero mean 
magnetic field is called a small-scale dynamo. The mechanisms of the small-scale 
dynamo action are different depending on magnetic Prandtl numbers Pm = Z//77, 
where v is the kinematic viscosity of the fluid and rj is the magnetic diffusion due 
to electrical conductivity of the fluid. For large magnetic Prandtl numbers, the self- 
excitation of magnetic fluctuations is caused by the random stretching of the magnetic 
field by the smooth velocity fluctuations (see, e.g., [151 HS1 HZl El UHl HH HH1 [201 EI]). 
This type of dynamo has been comprehensively studied in direct numerical simulations 
(DNS) of forced turbulence [221 1231 123 [25] and turbulent convection [23 [21]. The 
nature of small-scale dynamo for low magnetic Prandtl numbers is different, e.g., it 
is driven by the inertial-range velocity fluctuations at the resistive scale. The small- 
scale dynamo at low magnetic Prandtl numbers has been studied analytically (see, 
e.g., [2S1 EHl EQl Ell E2J E31 [31]) for a Gaussian white- noise velocity field (so called 
the Kazantsev-Kraichnan model) and numerically (see, e.g., [HI [351 EH]) in a number 
of publications. Since the magnetic energy is not conserved, the second moment of 
magnetic field has anomalous scalings [371 121] . 

The small-scale dynamo instability is excited when the magnetic Reynolds number, 
Rm, is larger than the critical magnetic Reynolds number, Rm cr . Analytical models 
based on the Kazantsev-Kraichnan model of a homogeneous, isotropic, non-helical 
and incompressible velocity field, yield Rm cr « 410 at very low magnetic Prandtl 
numbers [29] . Compressibility of fluid flow causes strong increase of the critical magnetic 
Reynolds number at Pm <C 1 (see [29J ) . Similar tendency also has been recently 
demonstrated in analytical study [21] at large Prandtl numbers. Direct numerical 
simulations of small-scale dynamo in [HI [35j [36] of the Navier-Stokes turbulence show 
that Rm cr is around 200 for small magnetic Prandtl numbers, and it is at three times 
larger than for the small-scale dynamo at large and moderate Prandtl numbers (see 
[22 | 1231 124"] ) . These DNS results at large and moderate Prandtl numbers are in agreement 
with different analytical models [101 HS1 1521 121] . 

The existence of the small-scale dynamo for a large number of turbulent spectra at 
large Prandtl numbers has been demonstrated in [21]. When Pm ~ 1 the small-scale 
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dynamo exists even in the regime of very large Mach numbers (see the DNS results 
in [25]). This study has also shown that, for low Mach numbers (~ 0.1), the ratio of 
the growth rate of turbulence driven by solenoidal and compressive forcing is about 30. 
However, for higher Mach numbers (~ 10), this ratio is about 2. 

The small-scale dynamo action is different from the turbulent induction effect that 
cases production of the anisotropic magnetic fluctuations by the tangling of the mean 
magnetic field by the velocity fluctuations (see, e.g., [381 EH EH SH 1121 03]). This effect 
cannot be described in terms of the small-scale dynamo instability. 

In spite of a number of studies of small-scale dynamo at low magnetic Prandtl 
numbers, there are some key questions that are subject of discussions in the literature. 
One of them is related to the scaling for the growth rate A of small-scale dynamo 
instability at low magnetic Prandtl numbers in the vicinity of the dynamo threshold. 
Our analysis performed in this study and even numerical solution of the dynamo 
equations for a Gaussian white-noise velocity field obtained in [H] imply that there 
are two different asymptotics for the dynamo instability growth rate: (i) in the vicinity 
of the threshold of the excitation of the small-scale dynamo instability and (ii) far from 
the threshold of the small-scale dynamo instability. 

Other issue studied here is related to an existence of the Golitsyn spectrum, k~ u / 3 , 
of magnetic fluctuations [38l [39] in the small-scale dynamo with low magnetic Prandtl 
numbers. This spectrum of magnetic fluctuations has been observed in the laboratory 
experiments [151 H6], in the large-eddy-simulations [17] and in the direct numerical 
simulation [HI [36] of the small magnetic Prandtl numbers magnetohydrodynamic 
(MHD) turbulence. In the present study we discuss conditions for the existence of 
the Golitsyn spectrum. 

The small-scale dynamo mechanism appears to be responsible for the random 
magnetic fields in the interstellar medium and in galaxy clusters [HI [131 113 EH [H]. 
A number of studies recently pointed out also the relevance of the small-scale dynamo 
to amplify small seed fields in galaxies and the intergalactic medium (see, e.g., 
[SH E21 [531 E3J EH] ) . In particular, DNS study in [51] demonstrated that in the presence 
of turbulence, weak seed magnetic fields are amplified by the small-scale dynamo during 
the formation of the first stars. Strong magnetic fields are generated during the birth 
of the first stars in the universe, potentially modifying the mass distribution of these 
stars and influencing the subsequent cosmic evolution (see [M])- It was also noted 
in [53] that the small-scale dynamo is very efficient during the formation of the first 
stars and galaxies. During gravitational collapse, turbulence is created from accretion 
shocks, which may act to amplify weak magnetic fields in the protostellar cloud. Such 
turbulence is sub-sonic in the first star-forming minihalos, and highly supersonic in the 
first galaxies. It was concluded in [53] that magnetic fields are significantly enhanced 
before the formation of a protostellar disk, where they may change the fragmentation 
properties of the gas and the accretion rate. 
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2. Governing equations 

Let us study magnetic fluctuations with a zero mean magnetic field at low magnetic 
Prandtl numbers. In sections 2-3 we use the Kazantsev-Kraichnan model [28] of the 
#-correlated-in-time random velocity field. Using this model allows to get the analytical 
results for the growth rate of the small-scale dynamo instability. The results remain valid 
also for the velocity field with a finite correlation time if the second-order correlation 
functions of the magnetic field vary slowly in comparison to the correlation time of the 
turbulent velocity field (see, e.g., pm [56]). The two-point instantaneous correlation 
function of the magnetic field can be presented in the form 

rW' 

(bi{t, x)bj{t, y)> = W(t, r)c% + — (5ij - r 4j ), (1) 

where W(t,r) = (b r (t, x) b r (t, y)) is the longitudinal correlation function, b r is the 
component of magnetic field b in the direction r = x— y, fy = TiTj/r 2 and W' = dW /dr. 
This form of the second moment (CQ) corresponds to the condition V • b = and 
an assumption of the homogeneous and isotropic magnetic fluctuations. Equation for 
the function W(r,t) derived in the framework of the Kazantsev-Kraichnan model of a 
homogeneous, isotropic, non-helical, incompressible and Gaussian white-noise velocity 
field, reads 

?*m = » w+Mw--%r, (2) 

at myr) m{r) 
(see [281129]), where 

1 2 2 r . . . 4 (\ 

myr) Km 6 mr V rn 

2/77 

«(r) =— f(r), f(r) = F(r)+rF'/3, 
r 

and Rm = uqIq/t] ^> 1 is the magnetic Reynolds number, uo is the characteristic 
turbulent velocity in the integral scale £q and F' = dF{r)/dr. Hereafter equations 
are written in dimensionless variables: length and velocity are measured in units of £o 
and uq- For a homogeneous, isotropic and non-helical (with zero mean helicity) and 
incompressible turbulent fluid velocity field the correlation function (r-Uj(x)-Uj(x + r)) 
is given by 

rF' 
Y 

The form of the continues function F(r) with different scalings in different ranges of 
scales is constructed using the following reasoning. The function F(r) = 1 — v^Re^ 2 is 
in the viscous range of scales, < r < £ v /£q, while the function F(r) = 1 — r 4 / 3 is in 
the inertial range of scales, £ u /£o < r < 1. At the boundary of these ranges, r = £ u /£o, 
these functions coincide, where £ v = £o/Re 3//4 is the viscous scale and £q is the integral 
scale of turbulence. 

The solution of Eq. ([2} can be obtained using an asymptotic analysis (see, e.g., 
[TU1 1291 [30]). This analysis is based on the separation of scales. In particular, the 



(TUi{x)uj{y)) = I 



(3) 
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solutions of Eq. (F2J) with a variable mass have different regions with different functions 
m(r), /i(r) and «(r). Solutions in these different regions and their derivatives can be 
matched at their boundaries. The results obtained by this asymptotic analysis are 
presented below. 



3. Asymptotic behaviour of the growth rate of magnetic fluctuations 

Let us discuss the asymptotic behaviour of the growth rate of magnetic fluctuations 
with a zero mean for small magnetic Prandtl numbers. We seek for a solution of 
Eq. ([2]) for the longitudinal correlation function of the magnetic field in the form: 
W(t,r) = exp(At) W{r). In the viscous range of scales, < r < £ u /£q, the function 
F(r) — 1 — ^/lier 2 and the equation for the function W(r) is given by: 

r W" + 4W' + ^-Pr m Re 3/2 rW = 0, (4) 
where W = dW{r)/dr. The solution of Eq. (J4]) is given by 



W(r) = r" 3 / 2 J 3/2 (y^^Re^rj » I - ^^ r 2 , (5) 

(see [29J ) , where J a {y) is the Bessel function of the first kind, and we have taken into 
account that W{r = 0) = I. 

In the inertial range of scales, £ u /£o < r < 1, the function F(r) = 1 — r 4//3 and the 
equation for the function W(r) is given by: 

f + - Rmr 4 / 3 ") W" + -(l + - Rmr 4 / 3 ) W 
3 / r V 9 / 

+ Rm(|r- 2 / 3 -^ W = 0, (6) 

where A is the growth rate of small-scale dynamo instability. In the range of scales, 
£u/£o < T £ v /£o the equation for the function W(r) reads: 

r 2 r + 4rW + — Rmr 4/3 W = 0, (7) 
27 

where £ v = £ /Rm 3 ^ is the resistive scale. The solution of Eq. ([7]) is given by 



W(r) = r~^J m lyn^ r V3 j w ! _ ^ r 4/a ; (g) 

(see [29J). On the other hand, in the range of scales, £ v /£q <^ r < 1 the equation for the 
function W(r) is given by: 

/ 27 A \ 

9 r 2 W" + 48 r W" + 52 r 2/3 W = 0. (9) 
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The solution of Eq. (jUJ) is 



^(^Cr-^JJ^r 1 / 3 !, (10) 
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(see [31]), where K a (y) is the real part of the modified Bessel function (Macdonald 
function) with a = («/2) \/39. This solution is chosen to be finite at large r, with 
positively defined spectrum, and it has the following asymptotics at scales r <C A~ 3//2 
(see ESI): 




W(r) = A 1 r~ 13/6 cos |y — In r + <p ] , ( I 1 

and at scales A~ 3 / 2 C r < 1 (see [32]): 



W{r) = A 2 r^ 7 / 3 exp I ^ r 1 / 3 j . (12) 

Here A\ and are the constants which are proportional to the constant C. 

In the range of scales r 3> 1, the turbulence is absent (F — > 0), 1/m = 2/3, 
«(r) = 4/mr and 

W(r) = A 3 r~ 2 (V\ + r" 1 ) exp(-Ar), (13) 



(see [29]), where v4 3 is a constant. 

The scaling for the growth rate of small-scale dynamo instability which is far from 
the threshold, is estimated as inverse resistive time scale: 

A-^-^Rm 1 / 2 , (14) 
t v to 

(see [2H]), where = (e ^) 1//3 is the characteristic turbulent velocity at the resistive 
scale, Wo = (^o) 1 ^ 3 an d s is the dissipation rate of turbulent kinetic energy. For the 
scaling A oc Rm 1 ' 2 , the condition A 1 ^ 2 r 1 / 3 ^> 1 implies r ^> Rm~ 3//4 . Note that the 
matching of the solutions (IE|), ( TL2l) and their derivatives at the boundary of their regions 
yields the dynamo growth rate (fT4"|) . 

1 /2 

However, the scaling, A oc Rm ' , is not valid in the vicinity of the threshold of 
the dynamo instability. Indeed, in the vicinity of the threshold when A — > 0, there is 
only one range of the solution of Eq. ([9]), i.e., A 1 ^ 2 r 1 / 3 <^ 1. In this range of scales the 
solution of Eq. is determined by Eq. ( TTT|) . The matching of the solutions OH]), CCD 
and fTl3|) and their derivatives at the boundary their regions yields the following growth 
rate of the small-scale dynamo instability in the vicinity of the threshold: 

A = /31nf^), (15) 



VRm cr / 

(see Appendix A), where /3 = 4/3 is the exponent of the scaling of the correlation 
function F(r) (i.e., it is the exponent of the turbulent diffusivity scaling). In Fig. 1 we 
plot the growth rate (1151) of small-scale dynamo instability versus ln(Rm/Rm cr ) in the 
vicinity of the threshold of small-scale dynamo instability. In the same figure we also 
show the numerical solution (squares) of the dynamo equation ([2]) performed in [44J (for 
the Kazantsev-Kraichnan model in inertial range of scales of fluid motions for Re ~ 10 s , 
and for Rm < 10 7 ), which demonstrates perfect agreement between the scaling ( Tl5|) . 
shown by solid line, and the numerical solution of the dynamo equation. 
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Figure 1. The growth rate of small-scale dynamo instability versus ln(Rm/Rm cr ) 
in the vicinity of the instability threshold: solid line corresponds to the scaling 
A = f3 In (Rm/Rm CI ) and squares are the results of the numerical solution of the 
dynamo equation for W(r) for the Kazantsev-Kraichnan model of velocity field with 
zero kinetic helicity in the inertial range of scales taken from Fig. 1 in [44] . 



Note that the solution of Eq. ((9]) determined by Eq. ( TlTl) . is generally a fast 
oscillating function at r < 1. However, for the first mode with the maximum growth 
rate, the spectrum of the eigenfunction is positively defined (see [ID])- Since this solution 
is only valid in the vicinity of the dynamo threshold, the second and higher modes are not 
excited. Therefore, the resulting spectrum of the eigenfunctions are always positively 
defined. 

In the present study we discuss only the regime of small magnetic Prandtl numbers. 
In the case of large magnetic Prandtl numbers and large fluid Reynolds numbers the 
dynamo growth rate far from the threshold is A ~ i?e 1//2 , i.e., it is determined by the 
Kolmogorov time scale (see, e.g., [39] EO] El])- On the other hand, for small magnetic 
Prandtl numbers the dynamo growth rate far from the threshold is determined by the 
resistive time scale. 

4. Magnetic fluctuations with the Golitsyn spectrum 

In this Section we study effect of magnetic fluctuations with the Golitsyn spectrum, 
A; -11 / 3 , [38| [39] on the small-scale dynamo with low magnetic Prandtl numbers. This 
spectrum can exist at the scales t v < r < £ v . Our goal is to determine the longitudinal 
correlation function W (r) that corresponds to the Golitsyn spectrum. To this end we use 
the induction equation for the instantaneous magnetic field H(t, x) in an incompressible 
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velocity field v(t, x): 

QTLT 

—— + (v • V)H = (H • V)v + 77AH. (16) 

We seek for the solution of Eq. ( fTB]) in the following form: 

H(t,x) = [B(t) + b(t,x)] exp(At/2), (17) 

where B(t) is the magnetic field in the scales which are much larger than the resistive 
scale £ v , while h(t, x) is the magnetic field in the scales which are smaller than £ v . We 
consider the magnetic dynamo regime, so that the total magnetic field H grows in time 
exponentially with the growth rate A/2. Now we average Eq. ([TBI) over the ensemble of 
fluctuations generated in the scales £ v <C £ <C £0, and subtract the obtained averaged 
equation from Eq. ([TBI) . This yields equation for the magnetic field b(t, x): 

— = (B • V)u + (rjA - A/2) b + b", (18) 

where v(t, x) = V(£) + u(£, x), V(t) is the velocity field in the scales which are much 
larger than the resistive scale £ v , while u(t, x) is the velocity field in the scales which 
are smaller than the resistive scale £ v , and b^ = V x (u x b — (u x b)) are the nonlinear 
terms. Equation (|T8|) is written in the frame moving with the velocity V(t). Using 
Eq. ( TT8]) and the momentum equation for the velocity u(t, x) we derive equations for 
the second moments of magnetic field hij(k) = (&i(k) bj(— k)) and the cross helicity 
tensor <7y(k) = (fe»(k) Uj(— k)): 

= - z(B-k) fo y (k) - ^(-k)] - (^ 2 + A/2) fc y + (19) 

^® = z(B-k) 4 (k) - (^ 2 + A/2) ^ + 4, (20) 

where ^(k) = («*(k) % (-k)), = (6f (k)6,(-k)> + (6,(k)6f (-k)) and = 
— i(k-B) /iy(k) + (6f (k)iij(- k)) + (6j(k)tt^(— k)). Here uf are the nonlinear terms 
in the momentum equation. Since we have already taken into account the exponential 
growth of the total magnetic field H, we can drop the time derivatives in Eqs. (fT9|) 
and (12"U|) . because the characteristic times of the variations of the correlation functions 
hij and ^ are much larger than the time A -1 . Since we describe magnetic fluctuations 
in the spatial scales which are smaller than the resistive scale £ rj , we may drop the 
nonlinear terms hfj and in Eqs. (fT9|) and ( 120]) in the case of large magnetic Reynolds 
numbers and low Prandtl numbers, because they are small in these scales. Therefore, 
the Eqs. (02$ and (J2D} yield: 

( V k 2 + A/2) 2 (6<(k) bj(-k)) = 2(B-k) 2 «k) Uj(-k)). (21) 

Now we introduce the normalized two-point correlation function w(r) of the magnetic 
field which is defined as follows: 

f(r) = (^n(x) H m (y)) = 1 + 7^ M*) My)>, ( 22 ) 

where w' = dw(r)/dr and r = |x — y|. We rewrite Eq. (12 ip in r space using the inverse 
Fourier transformation (i.e., we use the following transformation iki — > Vj). This yields 
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the following equation for the normalized two-point correlation function of the magnetic 
field: 



d 2 



2 d 



dr 2 ^ r dr 



w(r) 



'd 2 f(r) , 2df(r 



+ 



+ 3 a 4 , (23) 



3 y dr 2 r dr 

where a = (ARm/2) 1 / 2 , /(r) = f mm (r) and fij(r) is the two-point correlation function 
of the velocity fluctuations written in r space. Equation (1231) is written in dimensionless 
variables: length and velocity are measured in units of £q and uo, the growth rate of the 
magnetic fluctuations A is measured in units of uq/£q and the magnetic field is measured 
in units of B Q . We also take into account that ((B ■ V) 2 )/ = (1/3) (/" + 2/'/r). 

Solution of Eq. ( 123]) which satisfies the following boundary conditions: w(r = 0) = 3 
and w'(r = 0) = w"'(r = 0) = 0, reads 

3d 



w(r) = 3 — 



+ 



4 a 3 r 
55 Rm 2 



2(3 e 



6^ a ll/3 



[(ar — 1) exp(ar) + (ar + 1) exp(— ar)] 
(3 ar — 5) exp(ar) 7(2/3, ar) 



+ (3 ar + 5) exp(— ar) 7(2/3, — ar) 



(24) 



where w' = dw(r)/dr, C% is a free constant and 7(/3, x) = f3~ l x" exp(— x) M(l, 1 + x) 
is the incomplete gamma function which is related to the confluent hypergeometric 
function M(a,b,x). When ar <C 1, Eq. (1241) for the two-point correlation function w(r) 
is given by 



w (r) = 3 - d r 2 + — Rm 2 r 8/3 



1 + 
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238 



Rm Ar^ 



(25) 



The function w(r) is related to the longitudinal correlation function W(r), i.e., w(r) = 
3W(r)+rW'(r). Equation ( 126]) rewritten for the longitudinal correlation function W(r), 
reads 

RmRe 1/2 „ 1 „ s/! , f 9 
W(r) = 1 r 2 + — Rm 2 r 8/3 1 + — Rm A r z 
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(26) 



where £ u /£q < r < £ v /£q, the constant C\ ~ (5/3) RmRe 1//2 is determined by 
the matching of functions W(r) determined by Eqs. (jSJ) and (12^1) at the point r = 
£ v /£q. The scaling W(r) oc Rm 2 r 8 / 3 corresponds to the Golitsyn spectrum M(k) oc 
B%r]- 2 e 2/3 k- 11/3 [38], where e is the rate of dissipation of the turbulent kinetic energy 
and M(k) = (2/tt) / °° kr sin(kr) w(r) dr. It follows from the latter equation that the 
exponent q in the spectrum function M(k) oc k~ q and the exponent p in the scaling 
W{r) oc r p are related as follows q = p + 1. 

For small Prandtl numbers the constant C\ = RmRe 1//2 /3 is larger than (7 2 = 
Rm 2 /68, and since r<l, the first and second terms in the right hand side of Eq. (1261) 
dominate the behavior of W{r). This estimate implies that the third term in the right 
hand side of Eq. (126]) resembling the Golitsyn spectrum is negligible. This is the reason 
that the influence of the Golitsyn spectrum on the small-scale dynamo instability is 
minor. That is why it is so difficult to observe the Golitsyn spectrum in direct numerical 
simulations (DNS) for the small-scale dynamo with low magnetic Prandtl numbers. 
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Note that the existence of the Golitsyn spectrum of magnetic fluctuations requires 
a finite correlation time of the random velocity field, i.e., the solution for the small-scale 
dynamo with the Golitsyn spectrum does not exist in the framework of the Krachnan- 
Kazantsev model of the delta-correlated-in-time turbulent velocity field (see Appendix 
B). Indeed, for the derivation of Eq. ( )23l) we did not used assumption about the delta- 
correlated-in-time turbulent velocity field. One of the indications of the finite correlation 
time of random velocity field is already seen in Eq. (123]) . where the high-order spatial 
derivatives arise. The Kazantsev-Kraichnan model yields the dynamo equation with 
spatial derivatives not higher than the second-order spatial derivatives. On the other 
hand, it is well-known that even small yet finite correlation time of random velocity field 
causes the appearance of the higher-order spatial derivatives in the dynamo equations 
(see, e.g., [El EE!). 

The requirement of the finite correlation time of random velocity field for the correct 
description of the tangling magnetic fluctuations which have the Golitsyn spectrum even 
follows from the dimensional arguments. Indeed, the main balance in the induction 
equation for the magnetic fluctuations (B ■ V)u ~ DAh which yields the Golitsyn 
spectrum, can be rewritten in the following form: (b 2 ) ~ t 2 (£) [(u 2 )] 2 B 2 /D 2 . The 
latter equation implies the requirement of the finite correlation time of random velocity 
field for the correct description of the tangling magnetic fluctuations. The similar 
arguments are also valid for the k~ x spectrum of magnetic fluctuations generated by 
the tangling mechanism at low magnetic Prandtl numbers in the scales £ v <C £ <C £q 
(see 0DHHII12103!). 

We stress again that both magnetic fields, B(£) and b(t,x), are small-scale fields 
(in scales which are much less than the integral scale £q of turbulence). In particular, 
B(t) is the magnetic field in the scales £ v <C £ <C £ , while b(t, x) is the magnetic 
field in the scales l v -C £ <§C £ v - These fields belong to the same mode generated by 
the same small-scale dynamo mechanism. In this section we used two magnetic fields, 
B(t) and b(t, x), to describe interaction of the magnetic fields of different scales by 
the tangling of the field B(i) of the velocity fluctuations which produces additional 
anisotropic magnetic fluctuations with the Golitsyn spectrum. The latter mechanism is 
the turbulent magnetic induction that is different from the small-scale dynamo. 

5. Conclusions 

In this study we investigated some key issues of small-scale dynamos in random velocity 
field with large fluid Reynolds numbers, a zero mean magnetic field and low magnetic 
Prandtl numbers. Contrary to the claim in [44], there are two different asymptotics for 
the dynamo growth rate: in the vicinity of the threshold of the excitation of the dynamo 
instability [A oc ln(Rm/Rm cr )] and far from the dynamo threshold (A oc Rm 1//2 ). The 
influence of the Golitsyn spectrum on the small-scale dynamo instability is minor, and 
this spectrum of magnetic fluctuations requires a finite correlation time of the random 
velocity field. 



11 



We are grateful to Alexander Schekochihin for stimulating discussions who initiated 
our study of small-scale dynamo with the Golitsyn spectrum. We acknowledge the 
NORDITA dynamo program of 2011 for providing a stimulating scientific atmosphere. 



Appendix A. Growth rate of small-scale dynamo instability in the vicinity 
of the dynamo threshold 

Let us obtain the scaling for the growth rate of the small-scale dynamo instability in 
the vicinity of the dynamo threshold. In the range of scales, 4/4 < r < 4/4, the 
function rW'/W is given by 

^ = -iiw/3, ( A ,) 

[see Eq. OH])], while in the range of scales, 4/4 <C r < 1 the function rW /W is 

rW 13 /l3 / /l3 , \ , . v 

= \ — tan W — kir + (p , (A.2) 

W 6 V 12 IV 12 T y J 

[see Eq. fllip ]. On the other hand, in the range of scales r ^> 1 the function rW'/W is 
rW 1 

-2-VXr + (A.3) 



W 1 + v/Ar 

[see Eq. f lTB"]) ]. where we have taken into account that in the vicinity of the dynamo 
threshold A — > and A^ 2 r 1. 

Matching of the functions rW /W determined by Eqs. ( 1A.1I) and (1A.2I) at the point 
r = 4/4 yields the following equation: 



V39 1T1 31 /A . 

tan —luRm -„„ =-=. (A.4) 



This equation determines the function (po(Rm). Now we define the function = 
ifo(Rm = Rm cr ), where Rm cr is the threshold for the excitation of the magnetic 
fluctuations. It follows from this equation that 



v 39 , / Rm \ ... 

Matching of the functions rW'/W determined by Eqs. (IA.2j) and (IA.3j) at the point 
r = 1 yields 

fl2 /5 \ /A , 

tan ip = \ — - + A . (A. 6) 
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It follows from this equation that 



32 

A =37S M,) - (A ' 7) 



where we have also taken into account that in the vicinity of the dynamo threshold 
A —7- 0. Combining Eqs. (IA.5[) and (IA.7j) . we obtain the following scalinf for the 
growth rate of small-scale dynamo instability in the vicinity of the dynamo threshold: 
A = (4/3)ln(Rm/Rm cr ). 
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Appendix B. Tangling magnetic fluctuations in the delta-correlated-in-time 
velocity field 

The technique of path integrals for the delta-correlated-in-time velocity field allows us 
to derive the equation for the second-order correlation function, hij = (bi(t,x)bj(t,y)): 

= [L ik (x.)5 js + L js (y)5 ik + M ijks ]h ks + (B.l) 

(see for details [22]), where the turbulent component of magnetic field is b(t, x) = 
H(t,x)-B(t), 

U = -(l + — (B.2) 
3 3 V Rm/ lJ dx p dx p ' v ; 

1 M - A A f ° 2 A 9fm i 9 A 9 -4- 92 k m ^ 
7, M ijks — Oi k Oj s J mn — ik — — 0j s — — h — — , [D.6 

2 ox rn dy n dy s dx m dx k dy n dx k dy s 

d 2 f- 

Iij = BkBs dxi)f s > (B - 4) 

and f mn = (rw m (x)u n (y)). Multiplying Eq. ( IB. II) by r^rj/r 2 we arrive at the equation 
for the correlation function W(r,t): 

8 J^ = J-W" + Kr)W'-^\w + I, (B.5) 
at m(r) m(r) 

where I = B 2 (F" + 4F'/r) /3. In the inertial range the source term is / = 
52jB 2 r -2/3/27. This source term yields the following scaling of the correlation function 
W(r) oc r 4 / 3 in the range of scales, £ u /£o < r £ v /£o ■ This implies that the scaling 
of magnetic fluctuations caused by the tangling of large-scale magnetic field by the 
delta-correlated-in-time velocity field coincides with the scaling of turbulent magnetic 
diffusion F(r) oc r 4//3 . In Fourier space this corresponds to the k~ 7 ^ 3 spectrum of the 
magnetic fluctuations. This implies that the Golitsyn spectrum, k~ u ^ 3 , of magnetic 
fluctuations cannot be described in terms of the delta-correlated-in-time velocity field. 
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